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T— I ■ Abstract 
(N 

' Recent microscopic and phenomenological calculations of giant dipole resonances 
, for A < 56 nuclei are presented. The derived photodisintegration cross sections 
[ are exhaustively compared to the photonuclear data available to date. An accurate 
Qh' description of the data is found. Our new calculations are also compared with the 
Q . previous and widely- used estimates of Puget, Stecker and Bredekamp. The present 
calculations also include all the possible paths down the nuclear chart. The impact 
^ on the photodisintegration of ultra-high-energy cosmic rays (UHECR) is illustrated 
for a Fe source with typical energies of 10'^''~^^eV. At energies around 10^'' eV, 
the new cross sections are found to modify the UHECR photodisintegration rates. 
^ ' At energies around lO^^eV, it is recommended to solve a full reaction network to 
estimate the photodisintegration rate of the UHECR. 
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1 Introduction 



Cosmic-rays have been observed up to energies of ~ 3 10^° eV [1,2], which raises a double 
problem. First, one has to find an accelerator capable of reaching such ultra high energies. 
Second, it has been known for almost four decades that the interaction of the cosmic mi- 
crowave background radiation (CMB) with ultra-high-energy cosmic rays (UHECR) should 
lead to a sharp decrease in their flux above energies around 10^0 eV [3,4]. In the case of 
protons, which has given rise to the most detailed investigations, this cut-off is due to 
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pion photoproduction from CMB photons with energies above the reaction threshold in 
the proton rest frame. It has been proposed, however, that heavier nuclei could contribute 
significantly to the high energy part of the spectrum, as originally investigated by Puget, 
Stecker and Bredekamp (PSB) [5,6]. The propagation of nuclei in the intergalactic medium 
is also influenced by interactions with the background radiation fields, through photodis- 
integration reactions. 

In the nucleus rest frame, at typical UHECR energies of 10^^-10^^ cV, the CMB photons are 
boosted to energies in the range between a few hundreds of kcV to a few hundreds of McV. 
The interaction process between the UHECRs and the CMB is dominated by the giant 
dipole resonance (GDR) at photon energies below 30-50 MeV, and to a lesser extent by the 
quasideuteron emission for intermediate energies (between 50 MeV and 150 MeV) and the 
pion photoproductions at energies above 150 MeV [5,7]. In the original PSB model, two 
major approximations arc performed to estimate the intergalactic UHECR propagation. 
The first one concerns the total photoabsorption cross section which is parameterized as a 
simple Gaussian function [5,6], abruptly cut below the theoretical reaction threshold. The 
second one is based on the use of a reduced reaction network, involving only one nucleus for 
each value of A, to estimate the time evolution of the UHECR composition. More precisely, 
assuming that the /3-decay of the unstable nuclei produced by photodisintegration is always 
faster than the corresponding photoemission rate, a unique nuclear path is followed from 
the initial ^^Fe to the final protons [5] . 

The aim of this work is to study the impact of these two approximations on the prop- 
erties of the UHECR photodisintegration, and to provide the scientific community with 
new accurate determination of the photoreaction rates. Reaction model developments al- 
low now for an accurate and reliable systematic estimate of the photodisintegration rate. 
More specifically, photoreactions have been extensively studied in the field of nucleosyn- 
thesis, where phenomenological parameterizations of the photoabsorption cross sections 
have been optimized during the last decades, and large-scale microscopic predictions have 
also emerged [9,10]. New compilations of experimental photoabsorption data also help in 
defining the degree of accuracy of the present reaction models to predict the corresponding 
cross sections. Important progress has also been performed in the same field of nucleosyn- 
thesis to solve large reaction networks exactly in order to follow the time evolution of the 
composition of the material in given astrophysical sites. Similar tools can therefore be used 
in the field of UHECR to test the PSB approximation of a reduced network. 

Section 2 describes improved calculations of the photodisintegration cross sections, us- 
ing several phenomenological and microscopic descriptions of the dominant E'l-strength 
function. To test the accuracy of such models, they are systematically compared with 
the available photonuclear data for nuclei with A < 56, and with the PSB Gaussian pa- 
rameterization. Section 3 is devoted to the study of the impact of the newly-determined 
photodisintegration cross section on the propagation of a ^^Fe source. The approximation 
of the original PSB reduced network is also tested in Sect. 3 by comparing the UHECR 
photodisintegration obtained by solving the full reaction network exactly. 
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2 Calculation of the photodisintegration cross sections 



To describe the changes in abundance of the heavy nuclei as a result of the interaction of 
the UHECR with the CMB, a nuclear reaction network including all interactions of inter- 
est must be used. The chosen set of nuclear species are coupled by a system of differential 
equations corresponding to all the reactions affecting each nucleus, i.e. mainly photodisin- 
tegrations and /3-decays. The rate of change of the number density Nz,a of nucleus (Z, A) 
with charge number Z and mass number A can be written as 



— — — ^^Z+l,A^i3 + ^^Z-l,A^i3 



+-^^Z,A+l\^n + ■^^Z+l,A+l\^p + ■^^Z+2,A+4\^a 

\^U+2 , Ar .Z+2,A+2 , .Z+4,A+S 

-r^^Z,A+2'^'y,2n + ^^Z+2,A+2\,2p + ^^Z+4,A+8-^7,2a 

+Nz+l,A+2^j^p^~^^ + Nz+2,A+5^j^a^~^^ + Nz+3,A+5^j^a^~^^ 

-Nz,A [aI'^EA.^^^ 



(1) 



where X^'^ is the /3-decay rate of nucleus (Z, A) and Ay^ its photocrosion rate followed 
by the emission of a single neutron {x = n), proton (x = p) or a-particlc {x = a) or the 
emission of multiple particles such as 2n, 2p, 2a, np, . . . , including all open channels for a 
given photon energy distribution. 

The CMB photon density n(e) depends only on the UHECR Lorentz factor 7 = E/M(? 
(where E is the UHECR energy, and M its mass) [6] . The calculations of the CMB density 
as a function of the photon energy e in the nucleus rest frame show that photon energies 
overlap with the nuclear GDR for 7 ranging from 5.10^ to 10^^. In the nucleus rest frame, 
the photodisintegration rate X^^^ can be estimated from the cross section cr-y^j;(e) by 



A^,^ = j n{e) cr7,x(e) c de (2) 

where c is the light speed. 

All nuclei lighter than the seed nuclei and located between the valley of stability and the 
proton drip line must be included in the network. Under the most natural astrophysical 
assumptions, UHECRs are accelerated out of the ambient gas, possibly enriched in Fe 
close to neutron stars or depleted in metals (i.e. nuclei heavier than H) if significant pho- 
todisintegration occurs during the acceleration stage itself. Therefore, if nuclei are indeed 
present among the UHECRs, it is expected that they typically include the most abundant 
elements found in the interstellar medium, i.e. essentially lighter than Fe. The interac- 
tion of UHECRs with the CMB is thus expected to include all possible nuclei resulting 
from the photodisintegration of the heaviest species and therefore involve all stable and 
neutron-deficient unstable isotopes with A <^ 56. 
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One of the most complete compilation of photonuclear data is provided by the 2000 IAEA 
atlas (2000) [11]. Nevertheless, as far as elements of interest in the propagation of UHECR 
are concerned, only a limited set of photonuclear cross sections are known, namely the total 
photoabsorption cross section as a function of energy for about 10 nuclei and the integrated 
total photoabsorption cross section for no more than 16 nuclei. All the remaining rates must 
therefore be estimated on the basis of theoretical reaction models. 

2.1 The El-strength function 

The uncertainties involved in any cross section calculation are not so much related to the 
model of formation and de-excitation of the compound nucleus itself, than to the evaluation 
of the nuclear quantities necessary for the calculation of the transmission coefficients. The 
total photon transmission coefficient characterizing the probability to excite by photoab- 
sorption a compound nucleus excited state is obviously one of the key ingredients for the 
evaluation of the photoreaction rates. In the specific astrophysical conditions considered 
here, i.e for UHECR energies of 10^^~^^ eV, this function is dominated by the El transition 
which is classically estimated within the Lorentzian representation of the GDR. Experi- 
mental photoabsorption data confirm the simple semi- classical prediction of a Lorentzian 
shape at energies around the resonance energy Egdr- One the most widely used form of 
the £'l-strength function is described by the Brink- Axel Lorentzian model [12,13] 



where Egdr and Tgdr are the energy and width of the GDR, m is the nucleon mass and 
X — 0.2 is an exchange-force contribution to the dipole sum rule. 

The Lorentzian description is known to be less satisfactory at energies away from the GDR 
peak, and in particular fails to describe the low-energy experimental data, namely the 
radiation widths and gamma-ray spectra [14,16]. Various improvements have been brought 
to the Lorentzian form, mainly by including an energy-dependence of the GDR width 
capable of modifying the low-energy behavior of the ii^l-strength [8,14-16]. For this reason, 
the photon transmission coefficient is most frequently described in the framework of the 
phenomenological Kopecky-Uhl generalized Lorentzian model [16]. In this approximation, 
the GDR width of Eq. (3) is replaced by an energy-dependent width of the form r(£^) = 
TgdrI^^ + ^T^T"^]/ Eqdr) where T is the nuclear temperature and equals zero in the case of 
photoabsorption reactions. This model is the most widely used for practical applications, 
and more specifically when global predictions are requested for large sets of nuclei. It also 
requires the determination of the GDR peak energy and width to be predicted from some 
underlying model for each nucleus. For practical applications, these properties are either 
taken directly from experimental compilations (e.g [11,17]) whenever available, or obtained 
from a droplet- type model [18] or some experimental systematics [17]. 

The phenomenological Lorentzian approach suffers, however, from shortcomings of various 




(3) 
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sorts, and most particularly its lack of reliability when dealing with exotic nuclei or energies 
away from the GDR peak. For this reason, models of the microscopic type have been 
developed which are hoped to provide a reasonable reliability and predictive power for 
the E'l-strength function. Attempts in this direction have been conducted within models 
like the thermodynamic pole approach [17], the theory of finite Fermi systems or the 
Quasiparticle Random Phase Approximation (QRPA) [19]. The spherical QRPA model 
making use of a realistic Skyrme interaction has even been used recently for the large-scale 
derivation of the £^l-strength function. In such models, mean field calculations such as 
Hartree-Fock BCS (HFBCS)[20] and Hartree-Fock Bogoliubov (HFB) [21] are performed 
in order to describe the nucleus ground state. On top of these calculations, the QRPA is 
used to describe in a microscopic way the GDR. The linear response theory allows to predict 
the dipole strength of the excited nucleus. The only input of these mean-field models is 
the nucleon-nucleon interaction, including the Skyrme part, and the pairing component. 
Details can be found in [22]. 

Global HFBCS+QRPA and HFB+QRPA calculations for practical applications were de- 
veloped in [9,10] and shown to predict the location of the GDR in close agreement with 
experimental data and also to reproduce satisfactorily the average resonance capture data 
at low energies. These aforementioned QRPA calculations have been performed for all the 
8 < Z < 110 nuclei lying between the two drip lines. Note that the latest HFB-I-QRPA 
calculation is using the BSk7 Skyrme nucleon-nucleon interaction, initially derived to re- 
produce at best the measured masses [23]. 

2.2 The photodisintegration cross sections 

The photoreaction cross sections are estimated with the code named Talys [24] which takes 
into account all types of direct, pre-equilibrium and compound mechanisms to estimate the 
total reaction probability as well as the competition between the various open channels. The 
photoreaction cross section is estimated at energies up to 50 MeV. The calculation includes 
single particle (nucleons and alpha) as well as multi-particle emissions. All the experimen- 
tal information on nuclear masses, deformation and low-lying states spectra is considered, 
whenever available. If not, global nuclear level formulae, and nucleon and alpha-particle 
optical model potentials are considered to estimate the particle transmission coefficients 
and the nuclear level. Details on the codes and the nuclear physics input (ground state 
properties, miclear level densities, optical potential) can be found in the above mentioned 
references [25,24]. These various nuclear inputs are known to affect the photodisintegration 
rates mainly around the corresponding threshold energies, but have a relatively lower im- 
pact than the 7-ray strength function. To estimate the accuracy of the different approaches 
available for the evaluation of an £'l-strength function, the four models presented in the 
previous sub-section are considered, namely the Lorentzian [13], the generalized Lorentzian 
[16], the HFBCS+QRPA [9] and the HFB+QRPA [10]. In addition to the El component, 
the E2, Ml and M2 contributions are also included in the calculation of the 7-ray strength 
as prescribed by [17], but are not varied due to their lower impact in comparison with the 
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GDR. The quasideuteron process is also neglected in the present work due to the limited 
photon energy range [7]. 



2.3 Comparison with the photonuclear data 

The photorcaction cross sections estimated as explained above are now compared with 
available experimental data [11] for nuclei with A<56. It should be recalled that, even 
for stable nuclei, the data on such nuclei are scarce: a cross section measurement, such as 
(7,lnx) or (7,abs), is available at various energies for less than half of the stable nuclei 
with 12<A<56. Total photoabsorption cross sections around the GDR peak energy are 
available for 10 nuclei, while the integrated total photoabsorption cross section is known 
for 16 nuclei [11]. The exhaustive comparison performed to compare the predictions with 
the data is illustrated below. 

To test the model predictions, we first consider the (7,lnx) cross sections, i.c the photodis- 
intcgration leading to a single neutron emission, but possibly also to other extra cjcctiles. 
This set is of particular importance for our present study since the one-nucleon emission 
is the dominant process for photon energies in the 0-50 MeV range. The model predictions 
based on four different prescriptions for the £^l-strength function is compared in Figs. 1 and 
2 with experimental data for 8 nuclei, namely ^^C, ^^Na, ^°Si, ^^S, ^^Cl, ^^K, ^^V and ^^Mn. 
Globally, the calculations are found to describe the data accurately, from light to heavy 
nuclei, even though a few isotopes still suffer from imperfect description, such as ^^C. In 
most cases, however, the behaviour of the cross-section is remarkably well described, being 
confined within the experimental error bars over the whole energy range. For astrophysical 
applications, this degree of accuracy is quite satisfactory, especially in view of the remain- 
ing imcertainties pertaining to the origin of UHECRs (distribution of sources, injection 
spectra, initial composition) and to their propagation in extragalactic space, depending on 
the essentially unknown strength and topology of the magnetic effects. 

In some cases, the QRPA prediction is found to overpredict the data, especially at high 
energies. This is the case for the HFB+QRPA prediction of the ^^8(7, Inx) reaction, or for 
the HFBCS-I-QRPA results for ^°Si, ^^V and ^^Mn. On average, both phenomenological 
Lorentzian approaches reproduce better the experimental data. In the case of ^^C and 
^^Na, both the microscopic calculations and the Lorentzian-based predictions agree with 
each other. It should however be recalled here that the Lorentzian formulas make use of 
the experimental total photoabsorption peak cross section, peak energy and width for the 
cases presented here, while the QRPA models do not. 

In order to illustrate the accuracy of multi-nucleon photoemission predictions, the ^^V(7,2nx) 
measured cross section is compared to the four calculations in Fig. 3. The agreement is 
good and all models predict similar cross sections, except the HFB+QRPA model which 
overestimates the cross section. The errors on such measurement are relatively large and 
prevent one from considering this data as a strong constraint on the models. 
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Fig. 1. Measured photoabsorption cross sections (7,lnx) state, compared to the predictions of 

the four models: Lorentzian (dashed hnc), gcncrahzcd Lorentzian (soHd hne), microscopic HF- 
BCS+QRPA (dotted hne) and microscopic HFB+QRPA (dash-dot hne) 
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Fig. 2. Measured photoabsorption cross sections (7,lnx) state, compared to the predictions of 
the four models: Lorentzian (dashed line), generalized Lorentzian (solid line), microscopic HF- 
BCS+QRPA (dotted line) and microscopic HFB+QRPA (dash-dot line) 
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Fig. 3. Measured photoabsorption cross sections ^^V(7,2nx)compared to the predictions of the four 
models: Lorentzian (dashed hue), generahzed Lorentzian (sohd hue), microscopic HFBCS+QRPA 
(dotted Une) and microscopic HFB+QRPA (dash-dot hne) 

The (7,p) cross sections provide complementary insight on the accuracy of the models. The 
integrated cross section is available for about 15 nuclei of interest [11]. The description of 
the data is found to be reasonable, and not to differ from one model to another significantly. 
An average agreement within a factor of 2 is obtained on the integrated cross section. 

In view of the present comparison between the model predictions and the available ex- 
perimental data, the generalized Lorentzian model is adopted for further calculations of 

the UHECR propagation. Since the nuclei of interest for the propagation of UHECR are 
located close to the valley of stability (see Sect. 3.2), experimental data provide a relatively 
accurate systematic determination of the GDR properties and no significant deviation from 
the phenomenological parameterization of the Lorentzian type is expected. 

2.4 Comparison with the PSB parameterization 

The generalized Lorentzian cross sections are now compared with the PSB Gaussian pre- 
scription. This comparison is made for total photoabsorption data. Although (7,abs) cross 
sections are available for a limited number of nuclei (e.g. Ahrens et al. [26]), it allows a di- 
rect comparison of the 7-ray strength, regardless of the other nuclear inputs in the reaction 
calculation. The ^^Si and ^°Ca photoabsorption cross sections are compared in Fig. 4. 

In their original work [5], PSB calculate the total photoabsorption cross section, relying 
on the total integrated cross section, available for about one third of the nuclei of interest 
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Fig. 4. Measured total photoabsorption cross sections, compared to the predictions of the PSB 
models (solid line) and the generalized Lorentzian (dashed line) 

in the UHECR problem. This modelling of the GDR was performed up to an energy of 
30 MeV. Though not strongly different from the Lorentzian calculations, it is clear that 
the PSB agreement with experimental data is not as good as the one obtained with the 
Lorentzian parameterization. This may be explained by two features of the PSB modelling. 
First the PSB parameterization assumes a Gaussian shape for the GDR, which can differ 
significantly from the experimentally-confirmed Lorentzian shape. Second, the width of the 
GDR is frequently overestimated by the PSB parameterization, compared to the data. An 
exhaustive comparison between the PSB and the Lorentzian parameterization for the ~ 50 
nuclei used in the PSB model, shows also the general trend of a greater width (typically 2 
MeV) for the PSB parameterization. The large PSB width can lead to significant discrep- 
ancies, especially close to the one-neutron separation energy threshold located around 8 
MeV. The GDR contribution close to the neutron threshold plays an important role since 
it triggers the photodisintegration of UHECRs having just enough energy to be affected by 
the CMB photons while propagating throughout the universe. An accurate description of 
this threshold is thus necessary to derive the exact shape of the resulting turnover in the 
UHECR propagated spectrum. 

The larger width also leads to a larger integrated value of the cross section between and 
30 MeV than in the Lorentzian case. On the other hand the total integrated cross section 
between and 50 MeV is rather similar in both models since the Lorentzian prescription 
extends further up to 50 MeV while the PSB Gaussian prescription is limited to energies 
below 30 MeV. Note that both models predict a similar location of the peak energy since 
they consider experimental information or experimentally-based systematics. 



3 Impact on the UHECR photodisintegration 

The intergalactic UHECR photodisintegration is calculated on the basis of the rates (Eq. 2) 
derived from the Talys cross sections described in Sect. 2. All stable and neutron-deficient 
unstable nuclei with A<56 are included in the reaction network (Eq. 1). The photodisinte- 
gration of A<4 nuclei are not considered. Only the interaction with the CMB is considered. 
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Fig. 5. Evolution of the average mass number < ^4 > with respect to the distance of a ^^Fe 
source for three different Lorentz factors. Left: 7 = 710^ {E = 3.6 x lO^'^ eV); center: 7 = 2 lO^'^ 
{E = lO^ieV); right: 7 = 5 10^° {E = 2.61 x lO^^eV). 

the influence of the infrared background radiation being negligible in most cases of interest 
[6]. For illustrative purposes, we restrict ourselves to the case of a UHECR source made of 
^^Fe only. Other astrophysically relevant cases will be discussed elsewhere [27] . Note that 
a full propagation model for UHECRs (including pair production, statistical fluctuations, 
etc.) is not the purpose of this paper. Most recent calculations can be found, e.g., in [28,29]. 



3.1 Impact of the photodisintegration cross sections 



Fig. 5 shows the evolution of the average mass number < A > as a function of the distance 
from the ^^Fe source, calculated with the four GDR prescriptions described in Sect. 2. For 
a given source distance, < A > is the average value of the calculated nuclei abundances : 



where Ni is the number density given by Eq. (1). The full reaction network (Eq. 1) is solved 

at each time step, taking into account all the open photoemission channels, i.e (7,n), (7,p), 
(7, a), (7,2n), (7,2p), (7,2a), (7,np), (7,n«), (7,pQ;). In other words, the abundance of each 
type of nucleus is derived by taking into account the contribution of all the production 
channels, from the source nucleus downwards the table of nuclides, with the appropriated 
weight derived according to the corresponding cross sections. The values obtained are thus 
equilibrium values, representing the composition which would result from the propagation 
of an inflnite number of nuclei up to the time considered. For smaller UHECR samples, of 
course, Poissonian fluctuations of the numbers of each nuclear species are expected. 

Three illustrative values of the Lorentz factors 7 are considered. The value 7 = 710^ 

corresponds to an initial energy of the ^^Fe nucleus of E = 3.6 x 10^° eV. In this regime, 
only the lowest energy part of the i?l-strength overlaps with the photon density n(e). The 
photodisintegration rate is therefore sensitive to the position of the photoemission threshold 
located at the neutron separation energy {Sn ~ 8 MeV). The value 7 = 2 10^*^ corresponds 
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to a total energy of E = 10^^ cV. The photon density is maximum around the neutron 
separation energy, and the photodisintegration rate is sensitive to the location of the GDR 
peak energy. Finally, the large value of 7 = 5 10^° {E — 2.61 x lO^^eV) corresponds to a 
regime mainly sensitive to the energy-integrated photodisintegration cross section. 

Fig. 5 shows that iron nuclei with 7 = 7 10^ propagate up to distances 10 times larger than 
in the higher energy cases, before being stripped from all their nucleons. This is due to 
the small overlap between the photon density and the photodisintegration cross section in 
the I0W-7 case. In this regime, the distance of propagation is sensitive to the low-energy 
tail of the £^l-strength function and the use of different prescriptions leads to significant 
differences in the propagation distance. In contrast, the high-7 cases mainly depend on 
the GDR peak location or integrated photoabsorption, and for this reason the propagation 
distance is less sensitive to the photoreaction details. In addition, all calculations predict 
similar integrated photoabsorption cross sections since they all make use of the energy- 
weighted sum rule to normalize the dipole strength [30] . The spread observed between the 
different curves in Fig. 5 reflects the impact of the remaining nuclear uncertainties on the 
UHECR propagation. 

It should be noted that during propagation, the Lorentz factor of each fragment of the 
original nucleus remains essentially the same, since in a first approximation the total energy 
is evenly distributed among all nucleons. The whole propagation thus occurs at constant 
7 and the interaction regime (at the threshold, around the peak or through the whole 
range of the cross section) remains the same along the fragmentation process. It should 
be reminded, however, that additional energy losses occur in real astrophysical situations, 
due to e"'"e~ pair production as well as pion production. These important refinements are 
not included here because we focus on the specific infiuence of the nuclear cross sections, 
but they will be studied in detail in a forthcoming paper. 

The UHECR mass distributions are displayed in Fig. 6 and show interesting features. The 
standard mass deviation is defined as : 



All the i^l-strength parameterizations are seen to lead to similar predictions for the shape 
and magnitude of the standard deviation from the average mass number < A > as a 
function of the distance to the source. This overall mass distribution mainly depends on 
the relative competition between the various open photoemission channels, rather than on 
the absolute photoabsorption rate. This conclusion is valid irrespective of the energy of the 
source. At maximum spread, the average deviation from the mean value < A > reaches 
values as large as 7 mass units, which means that many diff'erent nuclei are found with 
comparable abundances as secondaries of the parent ^^Fe nucleus. In the specific case of 
7 = 7 10^, the mass distribution ranges from Ne to Fe isotopes at a distance D — 2.5 Mpc. 




(5) 
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Fig. 6. Evolution of the standard deviation of the A distribution with respect to the distance 
of a ^^Fe source for three different Lorentz factors. Left: 7 = 710^ {E = 3.6 x 10^° eV); center: 
7 = 2 10^° {E = 1.0 X lO^^eV); right: 7 = 5 10^° {E = 2.61 x lO^^eV). 

3.2 Photodisintegration path 



In Sect. 3.1, the UHECR propagation distance has been estimated making use of the 
full reaction network (Eq. 1). Previous calculations were based on the reduced PSB path 
illustrated in Fig. 7. In this approximation, only one stable isotope is considered per isobaric 
chain and the corresponding isobars (i.e. nuclei with the same A) are not affected by 
competitive channels. However, as shown in Fig. 7, about 85 nuclei are involved in the 
^^Fe photodisintegration at 7=2.10^° and numerous open channels including /9-decay can 
compete (the Lorentz dilation of time allows /9-unstable nuclei with half-lives of the order 
of the hour to survive over a Mpc scale, and thus have a chance to interact with a CMB 
photon). Most of the stable nuclei involved in the photodisintegration process have more 
neutrons than protons. Neutron emissions are therefore favoured and the corresponding 
unstable nuclei will /^^-decay towards the valley of stability. Note that we consider here 
that a given nucleus is involved in the reaction network if its calculated abundance amounts 
about 10% of the most produced one at any given time during the photodisintegration 
process. 

Significant differences can therefore be expected between our new calculation (Fig. 5) and 
the original PSB results based on the reduced path and the Gaussian parameterization of 
GDR strengths. In particular, as seen in Fig. 7, for heavy nuclei (A>45), about 70% of the 
nuclei are shortcut by the simphfied PSB path. For fight nuclei (A<45), less than 42% of 
the nuclei are bypassed. 

Fig. 8 displays the average mass numbers < A > with respect to the source distance for 
the three 7 regimes. For each regime, three calculations are shown, namely (i) the full 
network calculation based on the generalized Lorentzian rates (solid line), (ii) the reduced 
PSB path calculation with the same generalized Lorentzian rates (dotted fine), and (iii) 
the original PSB results based on the PSB path and Gaussian photodisintegration rates 
(dashed lines). The effect associated with the type of reaction network adopted is seen 
to be significant and to increase slightly with the UHECR energy. The simplified PSB 
path leads to a longer propagation distance with respect to the full network calculation. 
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Fig. 7. Nuclei involved in the photodisintegration process of ^^Fe nuclei for 7=2.10^°. Unstable 
nuclei are in shaded squares, and the PSB path is indicated by the arrows. The mass number of 
each nucleus is written in the corresponding square. 



independently of the nuclear input considered. This path effect is stronger for heavier than 
for hghter nuclei. For A < 45, the curves show similar slopes. For the heavy species, the 
major differences stem from the large number of nuclei excluded from the PSB path, while 
the full reaction network calculations show that many isobars contribute to the nuclear 
flow. Within an isobaric chain, the photodisintegration cross section is usually larger for 
high Z-values, so that nuclei on the PSB path propagate up to large distances. 



The differences between the dashed and dotted lines in Fig. 8 reflect the impact of the 
newly-determined photoreaction rates with respect to the widely used PSB rates. This 
comparison also confirms the previous conclusion that the cross section effect is attenuated 
at high energies due to similar integrated photoabsorption cross sections. However, the 
effect of the low-cncrgy i?l-strength around the threshold against particle emission remains 
significant, as seen in the low-7 case. Both the path and the cross section have an impact 
on the propagation distance, whereas at high 7 values, the path effect is the only one to 
remain. 
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Fig. 8. Evolution of the average A with respect to the distance of a ^^Fe source for three different 
Lorentz factor. Left : 7 = 710^ = 3.6 x lO^o eV); center 17 = 2 10^° {E = 1.0 x lO^^eV); right 
: 7 = 5 10^° {E = 2.61 x lO^ieV). 

4 Conclusions 



Four phenomcnological and microscopic models have been used to predict the photodisin- 
tegration rates for nuclei with A < 56 nuclei on the basis of the reaction code called Talys. 
The four models reproduce accurately the photonuclear data, namely the one-particle emis- 
sion, two-particle emission, and total photoabsorption data. The original PSB description 
of the total photoabsorption cross section often overestimates the GDR width, implying a 
larger GDR contribution at low energies with respect to experimental data. 

We studied in particular the case of a ^^Fe source and found that the new photoabsorption 
cross sections mainly modify the disintegration of UHECRs at relatively low Lorentz factors 
(~ 7 10^) , which is relevant to the shape of the cosmic-ray energy spectrum around the GZK 
suppression. The key ingredient proves to be a precise description of the GDR threshold. 
For higher Lorentz factors, details of the nuclear photoreaction rates do not impact much 
on the propagation distance. In contrast, it is found that solving a full reaction network 
can modify the results obtained with simpler nuclear paths, most particularly for the A 
> 45 nuclei. Large spreads in the mass distribution reaching ± 7 mass units are obtained 
when use is made of the full reaction network. 

The photodisintegration cross sections calculated with the generalized Lorentzian model 
are tabulated and made available to the scientific community via the nuclear astrophysics 
library at http://www-astro.ulb.ac.be. These include the (7,11), (7,p), (7,0;), (7,2n), (7,2p), 
(7,20;), (7,np), (7,nQ;), (7, pa) . . .for all nuclei with 12 < A < 56. 

The new approach to UHECR photodisintegration described in this paper will be applied to 
sensible astroparticle situations in forthcoming works, both using astrophysically motivated 
source compositions and taking into account the intergalactic magnetic fields. 
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